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Abstract 

In this paper, we prove a particle approximation, in the sense of the propagation of chaos, of a Lagrangian 
stochastic model submitted to specular boundary condition and satisfying the mean no-permeability condition. 
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1 Introduction 


In this paper, we construct a particle approximation of the following Lagrangian stochastic model (X, U) on a 
finite time interval [0, T], submitted to specular reflections at the boundary of a compact smooth domain 2? of R'*; 


Xt — Xq + f Ugds, Ut — Uo + f B[Xs; p(s)]ds + cWt + Kt, 
Jo Jo 

Kt =—2 ^ {Us- ■ n-viXs)) n-D{Xs)l^x^eav}, 


0<s<t 


, p{t) is the Lebesgue density of the law of {Xt, Ut) for t G (0, T]. 


( 1 . 1 ) 


The initial condition (Xq, Uq) is distributed according to a given probability measure po, and is independent to 
the R'^-Brownian motion {Wpt G [0,T]), nz> is the outward normal unit vector of the smooth boundary dV. 
We are considering Lagrangian stochastic model, this means that the dependencies in x of the coefficients in the 
velocity equation (1.1) are expressed as a conditional expectation with respect to the event {Xt = x}. Here the 
drift component B[x\ p{t)\ is a version of the conditional expectation E ]}){Ut) \ Xt = x]. Thus given a kernel b, 
B is defined for {x, j) G T) x L^{'D x R'^) as 


B[x-,j] 


Igd b{v)'Y{x,v)dv 
fjnd7{x,v)dv 


if 


/ j{x, v)dv ^ 0, 


0 otherwise. 


( 1 . 2 ) 


A particle approximation of Lagrangian Stochastic Models (SLM) like (1.1) has been studied without the confining 
jump term {Kt,t G [0, T]) in [4], and in the particular one-dimensional confinement case in [2]. The existence and 
uniqueness (in the weak sense) of a solution to (1.1) has been established in [3] for confined system in a smooth 
and compact domain V, with the help of PDE techniques. Furthermore the unique solution (X, U) satisfies the 
boundary condition (1.5), and the sequence of jump-times 


T„ = inf{T„_i < t < T; Xt G dV} for n > 1, tq = 0, 
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is well-defined and strictly increasing with n up to T, with the convention that inf 0 = T. 

Numerical algorithms for SLM are based on particle approximation methods (see e.g. [1] and the refer¬ 
ences therein). Here we give a first convergence result of a particle approximation of (1.1). We study the 
limit behavior of the interacting particle system = 1,..., N}, on a given proba¬ 

bility space (H, T, (J"*; t > 0), Q) endowed with independent copies {(Xg, Uq, (W^; t € [0,T])), z = 1,..., N} 
of (Xg, C/q, {Wt; t € [0, T])), defined as the solution to the following SDE system: 




= Ul,+ ! + aWl + K, 




i,e,N 
t 1 


(1.3) 


K 


i,e,N 


= -2 




i = l,...,iV. 


0<s<t 


where•= X is the empirical measure of the iV-particles, and is its marginal 

according to the two first canonical coordinates {xt,ut) at a given time t: ^ jji,e,Ny The 

drift [x; 7 ] is a smoothed version of B[x\ 7 ] in (1.2), with the help of a family of mollifiers (/>,(x) := e ), 
for some (j> € C^iV) such that </> > 0 and dx = 1. B^[x; 7 ] is defined for all x G 22 and all 7 in the set 

of probability measures on 22 x as 


(x,7) Se[x;7] 


fvxRd Hv)My)</>e(x - yhjdy, dv) 
fvxR<‘ lde{y)<t>e{.x - yh{dy, dv) + e 


(1.4) 


where Pe{y) = l{dist(y dv)>e} ctitoffs the support of 7 from a distance e to 922. 

The existence and uniqueness in law for the solution of (1.3) simply follow from Girsanov’s transformation 
and from the wellposedness of the confined Langevin process (i.e. the case 5 = 0, see Theorem 2.1 in [3]). This 
step only requires that 22 has a C^-boundary and that the support of /xg is included in 22 x 

Our main result is stated in Theorem 2.1: as the number of particle grows to infinity and the mollifiers param¬ 
eter e goes to 0, we prove that the particles (1.3) propagate the initial chaos with a limit law given by the solution 
to ( 1 . 1 ). 

Remark 1.1 (About the mean no-permeability condition on 922). In [3], we prove that the solution to (1.1) 
satisfies the so-called mean no-permeability condition: for x G 922, 


K[{Ufn-D{Xt))\Xt=x\ = Q. 


(1.5) 


Stochastic Lagrangian models have been introduced for complex simulation in Computational Fluid Dynamic 
( CFD). The mean-Dirichlet boundary condition (1.5) grounds the stochastic particle algorithm used to downscale 
simulations in CFD applications (we refer to [4][1 ] and the references therein for further details). 

Notice that the particle approximation of (1.5) (with the kernel b(u,x) = u ■ nxi{x)) in a neighborhood of 
922 is still an issue, that seems to require the continuity of the density of {Xt, Uf) over 22. Except in the one¬ 
dimensional case studied in [2 ], and to the best of our knowledge, such regularity result is unknown in the PDE 
literature on trace problems. 

Remark 1.2 (About the sequence of passage times on 922). When 22 = (0, -|-oo)xR‘^“^ and 5 = 0, the explicit 
expression of the joint law of (r„, ,n > 1) enables to control uniformly the confinement process ; t G 

[0, T]). For more general domain, we compensate the lack of such control by studying the trace problem for the 
density p. 

Notice that the estimate (3.6) in [2] on the upper-bound o/P(t„ < T) contains a mistake, claiming that this 
probability decreases with n uniformly in T. This shall be reformulated as follows: when the initial law po has 
its support in ((0, -Poo) xR"^”^) X R'^, there exists a constant C{T, a, fi*) depending on T, a and the distance 
fi* = sup{/3 > 0; supp{p,o) C ([/3, -boo) x R'*"^) x R"^} of the support of p,Q to ({0} x R'^"^) x R^*, such that, 
for all n > 3, 

P(r„<T)<C(T,a,/3*)^. (1.6) 

This clarification of the constant in front ofXj^P in the right hand side does not impact the results in [2], as we 
worked with fixed T. For completeness we give a short proof of (1.6) in Appendix A.2. 
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Notation. denotes the set ofpathsC([0, T]; 2?) xID)([0, xl])([0, T];R'^). Forallt G (0, T], we introduce 

the following sets: Qt := (0, t) xV x 

1]+ := \^{t,x,u) G (0,T) X dV x (it • nxi{x)) > 0} , 

I]“ := { {t, X, It) G (0, T) X &D X (it • ti-d^x)) < O} , 

:= { (t, X, It) G (0, T) X aP X (it • ni,(x)) = 0} . 

We set 12? I : = J-pdx. Denoting by dagu the surface measure on dV, the product measure onET:=S^USyUl]^ 
is dAsy := dt 0 dagT){x) 0 du. For a given positive weight function uj on R"^, we define the following weighted 
Sobolev spaces 

L^(w;Qt) := {/: OtR; ll/||i 2 (^.Q^) := f u}{u)f(t,x,u)dtdxdu<+oo}, 

J Qt 

yi(cc;Qr) := {/ £ C{[0,T]-L^tg-V x R'')); 

ll/lki0;QT) := ^max \\f{t)\\L2(u^.vxR'^) + II< +oo}, 

2.2(0;; E^) := {/ : ^ R; ||/|li 2 (^.s±) := fujiu)\{u-nvix))\fit,x,u)dX^^ < +oo}, 

l2(o;;2? X R'^) := {/: 2? X R'^R; ll/ll^s, := f w{u)f{x,u)dxdu<+oo}. 

’ JvxR<‘ 

A4 [E) denotes the set of probability measures on a measurable space E. When this is not ambiguous, we will use 
ll/llp for ||/||LP(i 5 ) with 1 < p < + 00 . 


2 Main results 

Hypotheses. From now on, we assume that the domain V, the distribution /to of (Xq, Uq), and the kernel b 
in ( 1 . 1 ) satisfy the following hypotheses ( 22 ). 

H-(i) dV is a compact sub-manifold of R'^. The initial measure /to has support in the interior of 2?xR'^ and 
/pxRrf l'a|2/io(da;, du) < +oo. /to has a density po in the weighted space L'^{uj] V x R'^) with 

u;(it) = (1 + |tt|2)T, for some a > (d + 3). 

H-(ii) b : —>■ R^ is a bounded continuous function and cr > 0. 

H-(iii) There exist Po’ Po ■ ^ in 2 .t(R+) , u i—>■ \/{l + |it|)Po(|it|) G L2(to;R‘^), and 

0 < PQ(|it|) < pq{x, It) < Po(|m|), a.e. on 2 ? x R'^. 

Notice that (22) are slightly more restrictive than the hypotheses in [3] for the existence: here b is assumed con¬ 
tinuous to simplify some weak convergence arguments, and the weight function ui is chosen in order to control 
Igd |u|3/w(it)dit. 

Theorem 2.1. Assume (22). Let P be the law on £ of {X, U, K) defined in (1.1), and let P*^W qJ^ 

{(X®’*^’^, 2T®’'^’'^), 1 < t < N} defined in (1.3). Then f’-chaotic; namely, for all {Ei, 1 < I < k}, 

k > 2, with Fi G Cb{C{[0, T]; 2?) x D([0, T]; R®^) x D([0, T]; R®^)), it holds that 


k 

lim lim (Pi 0 P2 ® ® Pfe ® 1 ® 1 ® P'^’'^) = TT(Pz, P). 

e—>-0+Af—)-+oo 


For the convenience of the reader, we summarize the results we need from [3], in the following. 
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Proposition 2.2 (see [3]). The solution {X,U) to (l.l) is unique in the subset of Ai{C{[0,T];'D) x I])([0, T]; 
that admits time-marginal densities {p{t); t G [0,T]) in L^{uj;'D x Moreover, (pit); t G [0,T]) solves in 
Vi{uj\Qt) the PDE 

{^2 

dtp{t, x,u) -Gu- Vxp{t, X, u) - — A„p(t, x,u) = - {B[x] p{t)] ■ Vup{t, x, u)) in Qt, 
p{0,x,u) = po{x,u) in VxW^, 

'Y~{p){t,x,u) =j^{p){t,x,u-2{u-nT>ix))nT>{x)) in 


where 7~'’(p) and 7 {p) are the trace functions of p, defined in L^(a;;E^) and L^(a;;Ey) respectively, and 
satisfies the following energy estimate 

\\p{^)\\l'^{lj-,'Dx9P) ^ ^ l|po|lL 2 ((^.x)xR<i) (1 + C'exp(Ct)), (2.2) 

JO 

where C > Q depends only on d, a and ||&||oo- In addition, p and its traces 7^(p) admit the following Maxwellian 
bounds: 

exp(a_t) (Go-(t) * £ 0(1 ■ 1))'"” (w) < p{t,x,u) < exp(a+t) (G^(t) * Pod ■ 1))'"^ (u) (2.3) 

andfordXsx- G S^, 

exp(a_t) (G^(t) + Pod • |))''” (u) < j^(p)(t,x,u) < exp(a+t) (G^(t) * Pod • 1))'"^ (u) (2.4) 

where Galt) is the centered Gaussian density function with variance a^t, * stands for the convolution product, 
a±, v± are constants depending only on T, d, a and ||&||oo <Jnd are such that a_ < 0, a+ > 0 and v± > 0. 

Solution of Equation (2.1) and the related notion of trace functions stated in Proposition 2.2 are understood in 
the weak (distributional) sense and we refer to [3, Definition 1.1 and Theorem 3.3] for a detailed existence result 
formulation. For the sake of completeness, let us mention that the existence of 7+ (p) and 7“ (p) is directly related 
to the notion of trace problems for kinetic equations (see references in [3]), and their construction is granted, on 
Sy and Ey respectively, by a density argument related to the solution space Vi(w, Qt) and the smoothness of 
dT>. Proposition 2.2 is clearly also true when the drift B[x\ ■] is replaced by its smoothed version Pe[x; •], or by 
a linear and bounded drift V{t, x). In the following, we slightly extend Proposition 2.2 for a generic linear drift 
V{t,x). 

Corollary 2.3. Let V G P°°((0, T) x V) and assume (H). Any weak solution (Xt, Up, t G [0,T]) to 
I X, = Xo + f*U.ds, 

\ Ut = Uo + fJV(s,X,)ds + aWt-2j:o<s<t(Us--nj,(X,))nj,(X,)lix.e07y} 


admits time-marginal densities {p{t); t G [0, Tj) such that p{t) is in L‘^{u!, T) X W^) for all t G [0, T], The traces 
7 “'"(p) and 7 “ (p) are in LF'{uj', E^) and L‘^{uj', E^ ) respectively, and such that for all f\ G CclTff), £ G Cc(Ey), 


E 


E 


h{^n,Xr„,U,^~)t{r„<T} 

nGN 

f2{Tn,Xr„,Ur„)t[T„<T} 


.nGN 


= / {u-nv{x))fi{t,x,u)y^{p){t,x,u)d\^,^ 

JY,t 


{u ■ n-D{x))f2{t, X, u)y (p)(f, x, u)dX^,^. 


( 2 . 6 ) 


Notice that Corollary 2.3 allows to extend the uniqueness result in Proposition 2.2 for the solution to (1.1) on 
the setX(C([0,r];P) x ©([O, T]; K'^)). 

The proof of Corollary 2.3 is postponed in the appendix. The rest of the paper is devoted to the proof of the 
propagation of chaos result. 

Although we give a particle approximation of the confined Lagrangian model, we are not able to use such 
approximation to construct a solution under lighter hypotheses than (H). In particular, we still have a deep us e of 
the PDE analysis of the Fokker Planck equation. The main difficulty resides in the uniform integrability result of 
the density traces, that we are able to show only with the strong Ma xwellian bound tool. 
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3 Proof of Theorem 2.1 


Equipped with the Skorokhod topology, £ is a Polish space. We denote by {Bt]t G [0, T]) the filtration associated 
to the canonical process {x{t),u{t), G [0,T]) off. 

The proof consists in the study of the double limits, first as TV tends to oo, next as e tends to 0. Mainly, we will 
detail the two following steps: 


Proposition 3 . 1 . Assume (H) and fix e > 0. The SDE 

U^ds, Ut = Uo+ f B, (s)] ds + aWt + , 

Jo 

Kt = -2 ^ {Ul- ■ nr,{Xl)) np(X:)l|xjeap}, h^{t) = Law{X^, Ut). 

0<s<t 



(3.1) 


has a unique weak solution, and we denote by the law on E of {X^, Then {P'^W- TV > 1} is P*^- 

chaotic; namely for all k > 2 and all (E), 1 < / < A:) of functions in Cb{£), 


k 


lim 

AT—>-+CO 


(Fi 0 • • • Ffc 0 1 0 1 0 • • • 01, P^’^) = Y[{¥f Fi). 

1=1 


(3.2) 


Lemma 3 . 2 . For any converging subsequence o/jP*^; e > 0} (that we still denote by P*^), the sequence of time- 
marginals densities {p'^(f) = P*^ o (a:(f), u(f))“^; e > 0} converges in L^(Qt) and in L^{u!;Qt) to the time- 
marginals densities {pit) = Q o [Xt, of the solution to (1.1), when e tends to 0. 

Proposition 3.1 has its analog in [2]. But now, the fact that the jump term is a finite variation process is not for 
free in the proof, as it is no more an increasing process. 

Notice also that even if we consider constant diffusion process, the mild-equation tool that we strongly used 
in [4] and [2] is useless here, as the controls we have on the semigroup derivative of the Lagrangian process are 
only for the L^-norm. 


3.1 The limit as N tends to oo 

Proof of Proposition 3.1. The wellposedness of Equation (3.1) directly derives from Proposition 2.2 (replacing B 
by B^) combining with Corollary 2.3. We only have to prove the propagation of chaos result. 

The verification of the Aldous’s criterion for the tightness of the family (P*^W. > 0) is a straightforward 

adaptation ofLemma4.4 in [2]. This ensures that the sequence = Law(^ ^ 

1} is tight on Al(f). 

We check that all limit points of TV > 1} have full measure on the set of probability measures un¬ 

der which the canonical process [x{f). u{t), k{t); t G [0, T]) satisfies (3.1). We denote by the limit of a 
converging subsequence of TV > l} that we still index by TV for simplicity. 

Eollowing Lemma 4.6 in [2], it is not difficult to see that, for 7r*^’°“-a.e. m G M.{£) with (m(T) := m o 
{x{t).u{t))~^-,t G [0,T]), the process 

Wt'.= — (^{f) — u{0) — k{f) — j i3£[x(s);m(s)] , f G (0,T], (3.3) 

is a R'^-Brownian motion under m. 

The remaining point is the identification of the jump process k that we detail in the following lemma. 

Lemma 3.3. The three following properties hold true 'K^'°°-a.e, m G Ai(£), m-a.s.: 

(a) For all jump times t G [0, T] of u, Au(t) = —2(u(t~) ■ nx>{x{t)))nx>{x{t)). 

(b) {k{t); t G [0, T]) is a finite variation process, and the related measure |fc| defined on [0, T] satisfies 

^{s > 0;x(s) G 922} ^ 

(c) The set {t G [0, T]; x(t) G 922} is at most countable. 
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The properties (b) and (c) above imply that {k{t)\t G [0,T]) is a pure jump process. In addition, since the 
paths of the process {u{t) — k{t)\ t G [0,T]) are continuous by (3.3), the jump times and the jump length of 
{k{t),t G [0,T]) and G [0,T]) are a.s. undistinguishable. Therefore (a) ensures that 7r*^’°“-a.e. m G 

M(£), m-a.s., k{t) = -2 • nx>(x(s)))nx.(x(s))l{ 2 ,,gai 5 } for all t G [0,T]. 

The uniqueness in law for the solution of (3.1) ensures that all converging subsequences of N > 1} 

tend to and enables us to conclude on the propagation of chaos property (3.2). □ 

Proof of Lemma 3.3. The proof mainly follows the proof of Lemma 4.8 in [2]. We only need to take care about 
points (Ij) as in the multi-dimensional case, the jump process k is no more an increasing process. 

For Pt = {p = {0 < ti < ... < f; < t}; / > 1}, the set of all partitions of the interval [0, f], the total variation 
process related to k is defined as 

i-i 

\k\{t) := sup ^ |fc(L-ri) - fc(L)l • 

First we prove that 7r'^’°°-a.e. m G A4{£), m-a.s., |A:|(r) < +cxd. We replicate some arguments of Sznitman [8] 
and introduce the sets 

Fm ■= |(a;, u, k) G £; |fc|(T) < M, J dist(a;(s), dT>) (i|fc|(s) = o|, := {m G M{£);m{FM) > 1 — r/} ■ 

Let us show that for all rj > 0, limM-i-+oo = 1- Since Fm is a closed subset of £ (see below), G^ is 

closed for the weak topology on M.{£). Therefore 

> limsup7r"’^(G2f) = limsup> 1 - p})- 

A^—>-+oo N—f-\-oo 

Denoting the complement of Fm on £, we have 

> 1 - 77 }) = 1 - > p}). 


Then applying two times the Chebyshev’s inequality, and using the exchangeability of the particles, 

> v}) < ■ 

Owing to Lebesgue’s monotone convergence theorem, we have 

z-i 


Eq [\k 


l,e,N 


It] = sup V E<i; 

m=0 




{xf‘’^edv} 




And, by the trace representation in Corollary 2.3, 

'-1 rtm+l 


Eq [|A:^’''’^|t] < 2sup ^ f [ \{u ■ nT)ix))\^-f {p^’''’^){s,x,u)dX^.^ 

'Pt ^_Q-/iTn Jd'DxR^ 

Since is bounded in L^(a;; uniformly w.r.t N, and — /R<i ^ + 00 , 



|(m • n-D(x))p7 (p^’''’^)(s,x,u)dAEj, 


< 



|(u • nT)(x))|a;(it)|7 {p^’'^X)\‘2[s,x,u)dX-s:T. 


< + 00 . 


It follows that lim sup^_^_|_(^ Eq [|iT^’'^’'^|T] < + 00 . Letting 77 tends to 0, we conclude that for 7 r*^’°“-a.e. m, 
™{Um>o^m) = 1 which means that \k\{T) < +cx 3 a.s. 
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We prove now that Fm is closed. Let us consider a sequence {(;■” = € N} in that subset, 

converging to = {x,u,k) in £ according to the Skorokhod topology; namely (see e.g. [5, Theorem 1.14, 
Chapter 6]) there exists a sequence {A„, n G N} of continuous increasing functions on [0, T] such that for all n, 
A„(0) = 0, A„(T) = T, lim„_^+oo suptg[o_T] |A„(f) - f| = 0, and 

lim sup |^"(A„(f)) — C(f)| = 0, and for alH G [0, T] lim |A^"(A„(f)) — A(^(f)| = 0. (3.4) 

n^+oo j-Q jij n —^+00 


As (3.4) implies that fc" o A„ and A fc” o A" converge respectively to k and Ak uniformly in [0, T], it further ensures 
that the sequence of measures o A„| converges weakly to d\k\ on [0,T]. Indeed, A" being a time change, 
let us first remark that the measure dj/c" o A„| coincides with the pushforward measure d(A“^tt|fc"|), such that 
foralls,f G [0,T], (A“itt|/c”|)([s, f]) = |fc”|([A„(s), A„(f)]). Since sup„ |fc"|(r) = sup„ |/c” o A„|(r) < M, 
there exists a converging subsequence o A„^|, f G N} with a limit denoted by 7 . From this sequence, let 
us further extract a subsequence o A„j;_ |, L G N} such that sup[Q o Xn^ — k\ < Then, for all 

continuous function / : [0, T] —>■ R, for the partition 0 = fo A A A ■ ■ ■ A < II = T of [0, T] such that 
\tm+i -tm\ <T/L,we have 




/ f(s)dlkl(s)- 

/ f{s)d\k'^^ 0 X^^\{s) 

Jo 

Jo 


< 


L-l 


+ 


+ 


/(s)d|fc|(s) - ^ f{tm) \k{tm+l) - k{tm)\ 

m—0 

-1 M-1 

\k{tm+l) - k{tm)\ “ f {tm) (A„j, (fm+l)) “ (A„j, (fm)) | 

m—0 m—0 

L-l 

^ f{tm) “ / f{s)d\k^^ O Xn^Hs) 

_n ^0 


L-l 


m—0 


The first and third terms in the right hand side tend to 0 as L —+00 by weak convergence. Since the second 
term is bounded from above by 2 max^g[g yj | /(f) |/L, we conclude that for any converging subsequence of |fc" o 
A„| we can extract a subsequence which converges to |fc|. This implies the weak convergence of \k'^ o A„| 
towards |fe|. Next, since t 1 —>■ r]n{t) := dist(x"'(f), 52?) and t p(f) := dist(x(f), 52?) are continuous and 
Jo Vn{s)d\k^\{s) = 0, 


V{s)d\k\{s) 


< 


[ V{s)d\k\{s)- [ p„(s)d|fc”|(s) 
'0 Jo 


? 7 (s)d|fc|(s) - / r]{s)d\k^ o Xn\{s) 


? 7 (s)d|fc” o A„|(s) - / ? 7 „(s)d|fc”|(s) 


The first term tends to 0, by the weak convergence of |fc" o A„| to |fc|. For the second one, using a change of 
variable in the second integral. 


1 

? 7 „(s)d|fc”|(s) = f 77 „(A„(A“ys)))d|fc”|(s) = f ? 7 „(A„(s))(iA“^tt|fc”|(s) = [ ? 7 „(A„(s))(i|fc”oA„|(s) 
Jo Jo Jo 


, we get 


? 7 (s)d|fc” o A„|(s) - / p„(s)d|fc”|(s) 


< sup |?7n(An(s)) - ?7(s)| M 
\s^[0,T] I 


from which we conclude that 


fo V(s)dlkl(s) 


= 0 , letting n tends to infinity. 


□ 


3.2 The limit as e tends to 0 

The tightness of {P*^; e > 0} can be shown again by replicating the verification of the Aldous’s criterion given in 
Lemma 4.4 of [2]. 


7 






















The main concern in that step is for the identification of the limit points. With Lemma 3.2, we easily check 
that any limit of a converging subsequence of e > 0} is a weak solution to the (1.1), as it satisfies the 
following martingale problem conditions 

(i) P° o (a;(0), 14(0), fc(0))“^ = po ® where Jq denotes the Dirac mass at 0 on R'^, by hypothesis. 

(ii) For all t G (0, T], P° o {x{t), u{t))~^ admits the positive Lebesgue density p{t). 

(iii) For all / £ the process 


f{x{t),u{t) - k{t)) - f{x{0),u{0)) - f (m(s) • Vx/(a;(s), u{s) - k{s))) ds 

Jo 

2 

{B [a:(s); p(s)] • Vuf{x{s),u{s) - k{s))) + y A„/(a:(s), m(s) - fc(s)) 
is a continuous P°-martingale w.r.t. the canonical filtration {Bt] t £ [0, T]). 

Indeed, as observed in [2], this is a direct consequence of the following convergence 



ds 


(3.5) 


lim limsup / \p\{x + h,u + 5) — p\{x,u)\dx du = Q, M t G (3.6) 

|/tM5K0 e_>.0+ JvxVA 

that can be immediately deduced from Lemma 3.2. 

(iv) P°-a.s., the set {t G [0, T]; x(t) G dV} is at most countable, and for all t G [0, T], 


k{t) = -2 ^ (u(s ) • ni5(a;(s)))iix)(x(s))l{a,(s)eax)}, 

0<s<i 

since we can reproduce all the arguments of Lemma 3.3., applying Corollary 2.3 again. 

Proof of Lemma 3.2. The existence of the time-marginal densities {p^{t)\ t G [0, T]) in L^(uj; Qt) follows im¬ 
mediately from Corollary 2.3. Adapting Proposition 2.2 to p'^ satisfies in Vi(w; Qt) the analog of the Fokker- 
Planck equation (2.1), replacing B by B^. Thus we observe that, for all e > Q, := p'^ — p satisfies in Vi (w; Qt) 

{ dtR^ + u ■ - y A,i?' = -V, • - B[--p]p) in Qt, 

I R'^{0,x,u) =0 inT> X 

[ X, u) = 'y~{R’^){t, x,u — 2(u ■ nT:{x))n'D{x)) in E^. 

Therefore, applying [3, Lemma 3.8] with B = B^[--,Pf], q{t,x,u) = "f~{R^){t,x,u — 2{u ■ nTi{x))nT>{x)), 
g = p^]-B[-,p])- V„p), for all t G (0, T], 

Jo 

= [ {B,[-,p^]p^-B[-,p]p)-{R^VuUJ + ujVuRn + ^^ [ (A,w)(i?^)^ 

J Qt Qt 


Since |A„w(m)| -f |Vua;(u)| < C{a, d)uj{u), using Young Inequality, it follows that 

< + (3.7) 

+ (^C{a, d) + ^ ||p(s) p'=(s)] - B[-, p(s)]) \\l2f^^.TxRf ^s- 
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Now, observe that, according to the Maxwellian bounds in Proposition 2.2, we can define finite positive constants 


M := sup ( / uj{v)\p(t,x,v)\^ dv] andm:= inf ( / p{t,x,v)dv], 

{t,x)e{0,T)xV yJs.'i / (t,x)G(0,T)xX> 


with m>0 and M < oo, from the Maxwellian bounds and H-(iii). Then 

rt 


[ ||p(s)(i3,[.;p^W]-i?[-;pW])|| 

Jo 


LHu-vxR'^)ds< M \B^[x;p'^{s)]-B[x;p{s)]\ dxds. 

J{0,t)xV 


By setting p(t,x) := p{t,x,u) du, bp{t,x) := J^ab{u)p{t,x,u) du, p'^{t,x) := p'^(t,x,u) du and 

bp^{t, x) := /jjd b{u)p'^{t, X, u) du and for + the convolution product, we have 


Beix-^p^it)] - B[x;p{t)] = 
= (j)^* {P^bp^){t,x) 


cj)e* (fdebp^) {t,x) bp(t,x) 


(j)e* iPePOit^x) + e p{t,x) 
p{t, x)-(j)e* iPeP^) {t,x)-e\ pe* Pebp^{t, x) - bp{t, x) 


{pe * iPePO {t,x) + e)p{t,x) 


p{t,x) 


pe*{PM){t,x) (P,*{p,p-P^p^){t,x)\ 1 , (Rijn jr\\(+ \\ 


p^ * iPepOit.x) + i 


p{t,x) 


Pe * {PebpO {t, x) f p{t, x) - pe* {PeP){t, x) - e\ pe * { Pebp){t, x) - bp{t, x) 


pe * {PepO{t,x) + I 


p{t,x) 


which gives 


l{ 0 ,t)xv 


\BPx] p'^{s)] — B[x; p{s)]\^ dxdsdu < 

m ^ 


{0,t)xVxR‘^ 
2 


p{t,x) 


Ip'^(s) - p(s)|^ dxds + Te 


foiTe := WPe * (/?eP) -Plli2((0,t)xi5) + “ ^Plli 2 ((o,t)xi5) + Comiugbackto (3.7) 


and 


using GronwalTs inequality, it follows that 




for C independent of e. Since p is in Lp{u}\ Qt) and lime_>.o /3e = 1 a.e. on 27, 7^ tends to 0 as e tends to 0. Hence 


liim^n IIi?' 


£^0 ll-rt \\Vi{uj-Qt) 


= 0 . 


□ 


A Appendix 

A.l Proof of Corollary 2.3 

From the Riesz Representation Theorem, it is sufficient to check that for alH G (0, T], there exists some a C > 0 
such that _ 

Eq yuj{Ut)P{Xt, Ut)] I < C||r/;||i 2 (i,xR<^), for all p G L\V x R‘^). (A.l) 

Without loss of generality, let us assume that p is nonnegative. Then 


^q[\/u}{Ut)p{Xt,Ut)] = / Eq ^Ju}{Ut)p{Xt,Ut)\{XQ,UP) = {x,u) pa[x,u)dxdu 
i'DxR'' 


< \B\ / Ec 
JVxR'^ 


y'ujiUPPiXt, Ut)\{Xo,Uo) = {x,u) -^^PoiluDdxdu 


= ||Po||i|27|Eq Zty/uj{vt)p{yt,vt) 


where {y, v) is defined as 


yt = Yo+ / Vsds, Ui = Vb + crlTt - 2 ^ • nr>(2/s))nr>(ys)l{y,Gax>}. 

Jo 


0<s<.t 
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with (yo,Vb) distributed according to and where Zt = exp(-i J*V{s,ys)dWs-^ f* \V{s,ys)\'^ds). 

The couple {yt,vt) is then distributed according to the density law h{t,u)dxdu := i^(Gcr(i) * 

deed, one can easily check that his a solution to the Fokker-Planck equation (2.1) in T^((0, T) x 22; 22^ (R‘^)) fl 
Vi{u!;Qt) with B = 0. The combination of Propositions 4.1 and 4.2 in [3] (in the case B = 0) gives Q o 
{yt, Vt)~^ = h{t, u)dxdu. Applying the Cauchy-Schwarz inequality, we get 


I5xR‘^ 


u)uj(u)h{t, u)dx du. 


Eq [v4;(t7t)V>(X,,[/*)] < \V\exp{'-i%^}^Eii[uj{vt)^P^yt,Vt)] < 

Observing that w(tt) < 2^ {uj{u — v) + a;(t;)), for all u,v £ R'^, we get 
uj{u)h{t,u) < 2^ ( f ljj{u — v)Ga{t,u — v)Po{\v\)dv + I w(z;)Gcr(f, u — u)Po(|'f|)d'f ) < 2'2+^C'(a) 


with G(a) = II wGcr(f, Oil oo II Poll ijthat allows to conclude on (A.l). 

Now we prove the probabilistic interpretation of trace integrals in (2.6). We consider the unique solution p 
in Vi(a;; Qt) and 7 ^(p) in L‘^{lo\ E^), of the following weak Fokker-Planck equation starting from pq\ for all 
t e (0, T], 4> G C“([0, f] X P X R"^), we have 


/ {u-nT>{x))j~^{p){s,x,u)(p{s,x,u)dX’s:T+ {u-nTi{x))^ {p){s,x,u)4>{s,x,u) dX^^ 

•2S+ 

= — (j){t,x,u)p{t,x,u)dxdu+ / (j){0,x,u)po{x,u) dx du (A.2) 


JVxR'^ JVxS.'^ 

+ J ^ds4> + u ■ yxfj) + V ■ \7 u 4> + {s,x,u)p{s,x,u) ds dx du 

and such that 7 ^(p) satisfy the Maxwellian bounds (2.4) (see [3, Proposition 3.14]). Then it i s straightforward 
to check that p{t) is also the density of {Xt, Ut) using the identification by mild-equation used in [3]-Proposition 
4.2. Now applying Ito’s formula to E[(/)(T, Xt, Ut)], combining with (A.2), one has 


E 


^ [<l>{s,Xs,Us- -2(U,- ■nT{Xs))nT{Xs)) - <l>{s,Xs,Us-))t{x^^av} 

0<s<t 

/ {u-nTi{x))^^{p){s,x,u)(j){s,x,u)dXT.T + / (m • ^^(x)) 7“(p)(s, x, u)(^(s, x, m) 

Jt± Jt.z 


Using the density between (P x R“) and Cc(P x R“) and the surjectivity of the application (j)£ Cc(P x 
(j)\^ G Cc(Et), we conclude on (2.6). 


A.2 Proof of the upper bound (1.6) on P(r„ < T) 

We consider the Langevin process on the probability measure Py,«, endowed with a Brownian motion B, 

Xt=y+ / Usds, Ut = v + Bt, 

Jo 

with y ^ 0, and the sequence of passage times 


T„ = inf{rn_i < t < T; Xt = 0}, forn >1, tq = 0. 


From the expression of Py,„(r„ G dt, |ut„ | G dz) given in Theorem 3, in Lachal [6], we obtain that, for any 
n > 2, 


pT p+oo , pt p+oc 2(z'^ 

Arn <T) = g{t - s, y, v; 0, u) exp ( - 


nt /* + Cx:) 


2{z^ + u^) \ ds du 


^+oo 


/o ^0 
7sinh(7r7) 


[Jo (2cosh(^)) 


_ TA. 


0 ^0 
/4uz 

'V“ 


) d'j 


(*+oo 


7sinh(7r7) 
/o (2cosh(fl))’ 


■K, 


i-y 


S 

Muz 
V s 


c?7 


(A.3) 


dz dt 
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where 


, ^ , 2^3 ' 

g{t - s, y, v; 0 , u) = — -^ exp 


'( 


6 t/^ 6yv 2{u‘^ + v'^)' 


) ( (t -^s)2 (^y + “ ®)^)) 


(i-s )2 ^V(t-s)3 (t-s )2 (t-s) J y{t-s) 

and Kij{a) = exp {—a cosh(i)} cos{jt) dt is the modified Bessel function. We then work with the expres¬ 
sion in (A. 3), using the following tricky identity, successfully used for similar computation in Profeta [7] 

sinh( 7 r 7 ) = sinh( —) f4cosh^( —) — . 

Assuming now that n > 3, for the first integral in A.3, we have the decomposition 


r+°“ 7sinh(7r7) 1 _ sinh(^) 


Jo ( 2 cosh(^))"-i s 


1 /■+ 


2"-i 

Furthermore, for all fc > 2, a > 0, we have (see [7], page 22) 
'•+“ sinh(^) '•+°° 


I 7^*7 ( ) 

7o s 


s 'cosh (^)"“3 
4uz^ sinh(^) 


(A.4) 


cosh(^)”“3 


d'y. 


° / <.»inh(e)exp(-<,c«h(«)) ^ ’I'l) ‘‘O- 


Since | sin( 76 >) I — ^nd sinh > 0 , for a = this yields to 


r+oc 

I 7^*7 (®) 

Jo 


sinh(z) I ^ 2 ^ 4uz 
— ’ s 

sinh(^) 


d'y 


cosh(^)"“3 

Using this last estimate, it follows that 


< 


p + OO 

/ a sinh( 0 ) exp ( — a cos 

Jo 


-too o 

e-¥d7)d6» = -e-“. 

TT 


/*+oo 


7 sinh( 7 r 7 ) / Auz 

K ' 


Jo ( 2 cosh( 7 ))- 




f+°° 7 sinh(7r7) ^ 

lo ( 2 cosh( 7 ))"-i s 77 


1 8 , Auz^ 

< t;— 5- — exp-. 

2"“3 jY V g J 


Therefore, coming back to (A.3), as exp(— 7 ^ ) exp(—7^) dz = exp{—^) dz < yTri, 

have ^ 

^y,v{_Tn <T) ^g{t-s,y,v,0,u)ds^ du^ dt. 

Let us now bound the integral 
/* + 00 


^+oo 

/ g{t- s,y,v,0,u)du 
Jo 

V3 f Qy 

:exp - 




6yv 2 u^ N — s) 


< — , = exp 

2yj2'K{t — s)3 


{t - s)3 {t - s)2 [t- s)) 2v/2 

32/" 


AV^(^-s) ( 1 / , 32/ ,2) 


(- -‘V 1 

V 2(t-s)^)' 


Thus, 


2{t — s)3, 

8/7r3 7 /■* 1 V3 




0 Jo - s)3 

8/7r3 1 7 -® y/3 


""P I - 2 (^J 


ds dt 


< 

— On-3 


r 1 r^- V3 / 3^7 ,,, 


that ends the proof, taking /3* = sup{/3 > 0; supp(po) C [/3, + 00 ) x R}, with the constant 




ds. 
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